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ABSTRACT 

We measure the alignment of galaxy ellipticities in the local Universe over a range 
of scales using digitized photographic data from the SuperCOSMOS Sky Survey. We 
find for a magnitude cut of bj < 20.5, corresponding to a median galaxy redshift of 
z w 0.1, and 2 x 10 6 galaxies, that the galaxy ellipticities exhibit a non-zero correlation 
over a range of scales between 1 and 100 arcminutes. In particular, we measure the 
variance of mean galaxy ellipticities, cr 2 (9), in square angular cells on the sky as a 
function of cell size and find it lies in the range, 2 x 10~ 4 > cr 2 {0) >^ 1 x 10~ 5 for cell 
side lengths between 15 < 8 < 100 arcminutes. Considering the low median redshift 
of the galaxies in the sample and hence the relatively low effective cross-section for 
lensing of these galaxies by the large-scale structure of the Universe, we propose that 
we have detected an intrinsic alignment of galaxy ellipticities. We compare our results 
to recent analytical and numerical predictions made for the intrinsic galaxy alignment 
and find good agreement. We discuss the importance of these results for measuring 
cosmic shear from upcoming shallow surveys (e.g. Sloan Digital Sky Survey) and we 
outline how these measurements could possibly be used to constrain models of galaxy 
formation and/or measure the mass distribution in the local universe. 

Key words: cosmology: observations - gravitational lensing - large-scale structure, 
galaxies: formation 



1 INTRODUCTION 



Over the last decade, considerable interest has been directed 
towards the measurement and analysis of galaxy ellipticities 
as a means to estimate shear induced by gravitational lens- 
ing (see e.g. Bartelmann & Schneider 1999). On the scale 
of galaxy clusters, this is now a well-established method, 
with tools to invert the shear pattern and measure the dis- 
tribution of Dark Matter. Recently attention has moved to 
larger scales to measure the cosmological weak shear signal 
of lensing by large-scale structure and a number of groups 
have made consistent detections (Bacon, Refregier & Ellis 
2000; Kaiser, Wilson & Luppino 2000; van Waerbeke et al. 
2000; Wittman et al. 2000). Despite the remarkable success 
of the shear analysis, up until recently very little attention 
had been paid to the prospect of intrinsic alignments mim- 
icking the gravitational shear signal. The most likely ways 
for this to occur is during the tidally induced spin-up of 
galaxies (Hoyle 1949), where the angular momentum axes, 
and hence ellipticities, are aligned, or through the alignment 
of galaxy and halo shapes. 

Over the last year, this problem has been addressed 
by a number of groups using a combination of numerical 
(Heavens, Refregier & Heymans 2000, Croft & Metzler 2001) 



and analytic methods (Catelan, Kamionkowski & Blandford 
2001, Crittenden et al 2001). Although these results are in 
rough agreement, a complete understanding of alignments is 
less secure, with the main problems lying in understanding 
the coupling of the tidal and inertial tensors of dark matter 
haloes, and the alignment of galaxies and haloes. On the 
observational side Pen, Lee & Seljak (2000) have recently 
claimed a weak detection of spin-spin correlations in the 
Tully catalogue. In this letter we have measured the variance 
of ellipticities of galaxies at low-redshift, where the effect of 
intrinsic alignments is predicted to be orders of magnitude 
higher than a lensing effect. These observations may help to 
distinguish between models for alignments. 

This paper is laid out as follows. In Section 2 we present 
the observational material used in the analysis. In Section 3 
the main analysis methods are introduced, and our results 
are presented in Section 4, where we discuss the significance 
of our results. Our conclusions are presented in Section 5. 



2 OBSERVATIONAL MATERIAL 

The observations analysed in this paper are taken from the 
SuperCOSMOS Sky Survey program (Hambly et al. 2001a). 
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This digitised photographic sky survey consists of Schmidt 
photographic plates (6° x 6°; plate scale, 1 cm = 11.2') 
and covers the entire southern sky (894 individual Schmidt 
fields) in two colours, bj and R. The material used in this 
analysis consisted of 436 Schmidt plates corresponding to 
~ 10, 000 square degrees.n. 

We have constructed an object catalogue (including 
both stars and galaxies) for the 436 fields by pairing the 
scanned bj and R plates. Apart from a small number of 
large overlap regions near the h boundary in the survey 
which were used for internal consistency checks (see Section 
3.4), we created a 'seamless' catalogue from the overlapping 
plates using the scheme described in Hambly et al. (2001a). 
This scheme attempts to include the image with the best 
parameters (and exclude the others) when there is a choice 
to be made for the same image appearing on more than 
one plate. Pairing the bj and R plates has the advantage 
over single-colour catalogues in that spurious objects on 
one plate are eliminated. This is particularly important for 
galaxy studies from Schmidt photographs since these defects 
are broken up into many co-aligned 'galaxies' by the image 
analyser. Such a large source of contaminants could poten- 
tially ruin any shear analysis. 

Regions around bright stars and blended (ie. multiple) 
images have been excluded from the object catalogue to fur- 
ther eliminate spurious and/or poorly parameterised images. 
Further details concerning image parameterisation, classifi- 
cation and photometry are given in Hambly et al. (2001b). 
Image parameters included in the final object catalogue 
generated for this study consisted of (for both the bj 
and R bands) celestial co-ordinates, local plate co-ordinates, 
second-order moments (semi-major/minor axes and celes- 
tial position angle, see Section 3.1), bj and R magnitudes, 
(bj-R) colour, image classification flag and stellarness index. 
The external reliability of the image classification on the 
J plates is demonstrated in Hambly et al. (2001b) as > 92% 
reliable for bj < 20.5 with completeness at around ~ 97%. 
Photometric accuracy for galaxies is around 0.25 mag. The 
internal consistency and accuracy in image ellipticity pa- 
rameters is demonstrated later in Section 3.4. 



3 ANALYSIS METHODS 

3.1 Measuring the ellipticities 

In order to measure the intrinsic alignment of galaxies, we 
first divide the survey into square cells of angular side length, 
9. We define the mean ellipticity of galaxies within each cell, 
6i, as 



ei.i 
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ei,j, 



3=1 



e, = e M + e 2>i , (1) 

where N is the number of galaxies in the i th cell and ex,j and 
ea,j are the ellipticity components of the j galaxy which 
are defined with respect to the cell axes as 
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Here, a and b are the semi-major and semi-minor axes 
of the galaxy and <p is the orientation of the semi-major axis 
with respect to the cell axes. The weighting scheme used 
is Wi — H{Ii — Ith) where h is the measured intensity in 
pixel i, I t h is a threshold intensity corresponding to the sky 
background and H is the Heaviside step function. That is, 
the SuperCOSMOS analyser measures second moments for 
each object detected as (Stobie, 1986) 

Qxx = [ d 2 6w i l{9){6 x -6 x f/ f d 2 6w t I{6), 



d 2 ewj(e)(e x ~6 x )(6 y -6 y )/ / d 2 e Wi i(e), 



d 2 0wil{e)(6 y -e y f/ / d 2 6wj(6), (3) 



where 9 X and 6 y are the centroids of the object and are given 
by 



d 2 ee x wj(e)/ / d 2 8wj(8), 



dr00 y wil(6)/ / d 8wil{8) 



(4) 



In terms of second moments, the ellipticity components of 
equation (2) are 
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The weighting scheme used corresponds to measuring un- 
weighted second moments within an isophotal threshold for 
all the objects detected. We have varied this threshold limit 
and have found the resulting measurements to be largely in- 
sensitive to the isophotal threshold used (see Section 4.2). 
To quantify the degree of alignment of the cell ellipticities, 
we calculate the variance of a on a given scale across the 
entire survey. Since this statistic has been used to measure 
cosmic shear (e.g. Bacon et al., 2000, Kaiser et al., 2000), 
we can directly compare the contributions of intrinsic and 
extrinsic galaxy alignments. In order to obtain an accurate 
estimate of e, for a cell, a number of corrections to the raw 
image catalogue are required. We now describe the sources 
of error on the ellipticity measurements of galaxies and the 
corrections we have applied to the dataset. 

3.2 Correction for PSF anisotropy 

There are several possible sources of error which could po- 
tentially compromise any shear analysis. These distortions 
result in spurious ellipticities for all the objects detected and 
must be removed before a shear analysis can be performed. 
Firstly, there may be slight astrometric distortions present 
in the dataset due to emulsion shifts in the photographic 
plates. However, these should be negligible — the dataset 
we have used has very precise astrometry and has been used 
as a standard for making astrometric corrections in other 
shear analyses (e.g. Gray et al., 2001). We have therefore 
not corrected for astrometric distortions in this analysis. For 
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Figure 1. Ellipticity fields for raw (LHS) and corrected (RHS) 
stars for three of the bj-band SuperCOSMOS fields. Stars are 
binned into cells of a side 10' and smoothed with a Gaussian 
with a smoothing scale 15'. For each plot, the length of each 
vector drawn is 25 cm times the measured cell ellipticity. The 
average ellipticity in a cell in the raw fields is e si 10 -2 , while the 
average ellipticity in a corrected cell is e zz 10 -4 . These plots are 
typical for the PSF anisotropy distributions on a plate and the 
corresponding R-band plots are generally very similar. 
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Figure 2. Ellipticity distribution per cell (stars on the left, galax- 
ies on the right) for the same three bj-band fields shown in figure 
1. The dots are the raw stars/galaxies, the crosses the corrected 
values. Each cell has an ellipticity estimated from equation (1) 
and has an uncertainty given by equation fl!5h. Average stellar 



distortions in the raw fields of e\ ~ e2 



10" 



are corrected 



to produce flattened stellar ellipticity distributions with residual 
cllipticities, e\ ~ e~2 ~ 10 -4 . Similar "before" and "after" distri- 
butions are found for the R plates. 



a detailed discussion of the astrometric accuracy of the cat- 
alogue, see Hambly et al., 2001c. 

The first correction we have made to the dataset is 
a correction for anisotropy in the point spread function 
(PSF). The are several sources of error contributing to the 
PSF anisotropy. Atmospheric dispersion along with emul- 
sion shifts and mechanical plate distortions, are minor con- 
tributors to the PSF anisotropy. By far the most important 
factors, however, are tracking errors and field rotation (for 
an analysis of the alignment, pointing accuracy and field 
rotation of the UK Schmidt Telescope, see Wallace & Trit- 
ton, 1979; further details of the telescope optics are given in 
Wynne, 1981). These effects combine to produce a system- 
atic PSF anisotropy pattern across each plate which needs 
to be corrected for before the galaxy ellipticity measure- 
ments can be trusted. We have done this by comparing with 
the ellipticity field for the stars. Stars should have no in- 
trinsic ellipticity and so the measured stellar ellipticities are 
due to the PSF anisotropy. As shown in Kaiser, Squires and 
Broadhurst, 1995 (hereafter, KSB), the perturbation to the 
galaxies' ellipticity components is given by 



Se a 



yf3Pf3, 



(6) 



where P^ g is the "smear polarizability tensor" which de- 
scribes the response of the individual galaxy images to the 
PSF anisotropy, p a which in turn can be measured from the 
stars. The analysis described in KSB, which is for a gen- 
eral weighting of quadrupole moments, is to measure the 
quantities, P^p for each individual galaxy and p a from the 
foreground stars and to apply the correction in equation M 
to the individual galaxies. However, for the essentially un- 
weighted moments that we have made use of in our analysis 
(see Section 3.1), P^g is diagonal with Pfj = P22 being a 
measure of the inverse galaxy size and the KSB correction 
reduces to 



Se a 



ixx + • 



xx T" ' 



(7) 



where the Q*'s and the Q's are stellar and galaxy second 
moments respectively and e* a are the measured stellar ellip- 
ticity components. In terms of semi-major and semi-minor 
axes, this becomes 
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Figure 3. The linear correlation coefficients, C\ and C2 of the 
cell-averaged stellar and galaxy ellipticity components as de- 
scribed in the text, after correction for PSF anisotropy. The 
shaded region is the 3cr significance level. 




Figure 5. C\, C2 and the 3cr significance level after excluding all 

9 g ) from the catalogue. 



easterns) 



e 1 (stars) 
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Figure 4. The observed anti-correlation of corrected galaxy cell 
ellipticity components (ei on top, e2 on bottom) with the corre- 
sponding observed stellar ellipticities. 
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Figure 6. Cell-averaged galaxy ellipticity components (ei top, 
e2 bottom) plotted as a function of the corresponding stellar el- 
lipticities after applying the galaxy size cut described in the text. 
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a 2 + b 2 



(8) 



where again, the superscript refers to the stars. We have 
applied the correction in equation (H) to individual galaxy 
ellipticities using cell-averaged measurements for the "stellar 
size" (a* 2 + b* 2 ), and where the e* 's are the average stellar 
ellipticity components in the cell. The size of the cells used 
for this correction was 10 x 10 arcmins. 

To perform the correction, the plates are gridded up and 
an average stellar ellipticity is calculated for each cell in the 
grid. This ellipticity field is then Gaussian smoothed across 
the plate. Different smoothing scales were investigated in or- 
der to find an optimum value. The resulting measurements 
however did not change significantly for smoothing scales 
between 10 and 30 arcmins. The individual galaxy elliptici- 
ties are then corrected according to equation M. If the grid 
used for smoothing is too small, a coherent distortion pat- 
tern cannot be made and an essentially random pattern is 
generated. Figure 1 illustrates the bj band stellar elliptici- 
ties on three plates before and after this correction. A strong, 
coherent distortion of e ~ 10~ 2 is corrected to produce a ran- 



dom ellipticity pattern with mean e ps 10~ 4 . Figure 2 shows 
a scatter plot of the ellipticities (ei vs. e-i) for both stars 
(LHS) and galaxies (RHS) before (dots) and after (crosses) 
correction. Each point corresponds to a cell, with an ellip- 
ticity value given by equation (1), and an uncertainty given 
by equation (113) . 

To test the success of the PSF correction further we 
calculated the linear correlation function of the cell-averaged 
stellar and galaxy ellipticity components (c.f. Bacon et al., 
2000), 



C, 



(e* a e a ) - (e*)(e a ) 
cr(e£)cr(e a ) 



(9) 



where e* a and e a are the cell-averaged stellar and galaxy el- 
lipticity components respectively and <j(e* a ) and cr(e a ) are 
the errors on those two quantities as measured from the 
data. In Figure 0, we plot C\ and C2 as measured from the 
total data set after correction for PSF anisotropy along with 
the 3cr significance level. It can be seen from this figure that 
there is a clear anti-correlation between the corrected galax- 
ies and the raw stars and furthermore, it is a S> 3u effect, 
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especially at large bin sizes. This anti-correlation is clearly 
seen in Figure H where we plot the corrected galaxy and 
raw stellar ellipticity components for a bin size of 6 — 70 
arcmins. This anti-correlation is the same effect seen in the 
cosmic shear analysis of Bacon et al. (2000) and is due to 
an over-correction of the PSF anisotropy for small galaxies. 
To remove this effect, we have imposed a cut on galaxy size, 
only taking galaxies with 6 g > 0.8 (0 g ), where 8 g = a 2 + b 2 
is the individual galaxy size, and {0 g ) is the mean galaxy 
size over the whole survey. Only galaxies satisfying this con- 
dition have been used in the final analysis. The correlation 
functions, equation (9), calculated after applying this cut, 
are shown in Figure 5. It is clear from this figure that the 
star-galaxy anti-correlation is no longer significant at the 
3cr level and this can be seen by eye in Figure g| where the 
corrected galaxy and raw stellar ellipticity components are 
plotted after applying the cut on galaxy size, again for a bin 
size of 70 arcmins. 



3.3 Seeing correction 

After correcting for PSF anisotropy, we model all other 
sources of error (atmospheric turbulence, wind shake etc.) 
as random effects. Collectively called seeing, we assume also 
it is isotropic and we apply a stochastic correction to remove 
its effects. The effect of seeing on the galaxy ellipticities is 
to circularize the images, causing a decrease in ellipticities. 
The seeing across all of the plates is typically 2". This is 
comparable to the size of a galaxy near the magnitude limit 
of our survey (bj < 20.5) and thus, needs to be corrected 
for before the galaxy ellipticities can be trusted. 

We assume that the semi-major and semi-minor axes of 
each galaxy transform under the effect of seeing as 



/2 _ 2 , 2 



u' 2 7 2 , 2 

b J = °J + r i ' 



(10) 



and (a,j,bj) are the axes of the j galaxy 
before and after the effect of seeing respectively and rt is the 
average seeing in the i th plate. The average galaxy profile as 
measured using the SuperCOSMOS machine deviates from 
Gaussianity only very slightly (see Knox et al., 1998) and 
so equation (lid) should describe the effect of seeing on the 
galaxy images reasonably accurately. 

The effect of seeing on the measured ellipticities is 



JijS-o 



-a, j 



(11) 



a "seeing factor" which is, in terms of the observed semi- 
major and semi-minor axes, 

„2 



/« = 1 
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(12) 



Here, a and b are the axes of the galaxies after correction 
for PSF anisotropy and the uncertainty on individual galaxy 
sizes was too large after this correction to use equation (111]) 
directly. Instead we estimated an effective galaxy size for 
all plates for a given flux cut. We have done this by look- 
ing at plate overlap regions. We have applied the correc- 
tion in equation (FLjLI) to galaxies in these regions and have 
"fine tuned" the effective galaxy size so that consistency is 
achieved in the measurement of both individual and binned 
galaxy ellipticities between the overlapping plates (see Sec- 
tion 3.4). In the case of bj < 20.5, the optimum effective 
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Figure 7. Individual galaxy ellipticities measured from overlap- 
ping SuperCOSMOS fields. The left-hand panel shows the ellip- 
ticity measurements before the seeing correction is applied. The 
plate plotted on the horizontal scale has greater seeing than the 
plate plotted on the vertical axis. The right hand panel shows the 
ellipticities after the correction. 
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Figure 8. Binned cell ellipticities, ej measured from the same 
plates as shown in Figure 3. The left-hand and right hand panel 
are the measurements before and after the seeing correction is 
applied respectively. 



galaxy size was 2.5", which closely matched the directly 
measured mean galaxy size. This was then used to correct 
all the plates, using the measured seeing from each plate. 



3.4 Internal consistency tests 

We have tested the seeing correction of the previous section 
by comparing ellipticity measurements in the plate overlap 
regions. The ellipticity measurements for a galaxy in an over- 
lap region as measured on the two overlapping plates should 
agree to within the limits of the measuring process (the mea- 
surement errors are dominated by noise on the original pho- 
tographs) assuming that we have corrected for the effect of 
seeing accurately enough. 

Figures |7J and M show the correlation between differ- 
ent plates in the overlap region, after correction for PSF 
anisotropy, for individual galaxies (Fig. f7J) and binned cells 
(Fig. H). The left-hand-side of Figures [fl and M show the 
correlation before correcting for seeing. There is an appar- 
ent decrease in the ellipticities of the galaxies plotted on 
the horizontal axis, where the seeing on the plate is larger. 
The right-hand-side of Figures m and M show the correlation 
after correction for seeing. As we are only scaling elliptici- 
ties the scatter is slightly increased, but the correlation is 
significantly greater. These plots are typical for the overlap 
regions. 
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Figure 9. LHS: Scatter plot of the SuperCOSMOS galaxy el- 
lipticities on a plate against those from the APM galaxy survey. 
RHS: The same as the LHS but for binned cells. 



3.5 Comparison with APM Sky Catalogue data 

To check our results further, we have performed a compar- 
ison test for similar alignments in APM dataFI in order to 
test for any systematic effects introduced by the SuperCOS- 
MOS scanning procedure. We took a J-R paired star/galaxy 
catalogue from the APM survey for one field and paired it 
up with the SuperCOSMOS data. Since the APM data is 
derived from a copy of the plate that SuperCOSMOS has 
scanned, any discrepancies in the scanning procedure should 
be highlighted by this test. Although the measurement of 
galaxy shapes for the APM machine is different, using den- 
sity measurements rather than the intensity measurements 
of SuperCOSMOS, the magnitude cut is low enough that 
the difference should not matter. Figure |9| shows a scatter 
plot for individual galaxy ellipticities for a plate on both 
SuperCOSMOS and APM surveys, and for the binned data, 
indicating a strong agreement between the surveys. 

Having checked for internal and external consistency 
in the ellipticity catalogue, we now turn to estimating the 
variance of the ellipticities in cells of varying scale, across 
all of the plates. 



3.6 Estimator for the ellipticity variance 

The variance of the cell ellipticities can be expressed as the 
sum of the contributions from all possible sources of align- 
ment (e.g. Bacon et al., 2000) 
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"tot 



2 
"lens 



2,2 
"int ' "noise 



i - 

+ a s 



(13) 



where we have included the contributions from lensing, in- 
trinsic alignments, shot noise and systematics. We assume in 
what follows that the lensing signal is negligible compared to 
that from intrinsic alignments for the median redshift of the 
galaxies in our sample (z fa 0.1). The noise term, a noise is 
due to intrinsic scatter in galaxy ellipticities and the random 
error in the measurement of the galaxy ellipticities. Since we 
are averaging over a very large area of sky (436 plates cov- 
ering ps 10,000 sq. degrees with fa 1.7 x 10 6 galaxies), we 
can expect to beat this term down statistically. The final 
term in equation ([13) is due to systematic sources of er- 
ror. In Section 3.4 we demonstrated the internal consistency 
of the catalogue, indicating that contributions to this term 



from variations in th e ca talogue are small. Further tests are 
presented in Section 4.2 



We have used a minimum variance estimator for the 
intrinsic ellipticity variance, a 2 nt due to intrinsic alignments 
in excess of the noise, 



X>i(e?-M) 



(14) 



where Mi is the random noise on the estimated cell ellipticity, 
d, in the i th cell, and Wi is an arbitrary weighting factor. 
We assume that the mean cell ellipticity components are zero 
(see Section 4.2). For the moment we will ignore systematic 
terms. Following the analysis of van Waerbeke et al. (2000), 
we estimate the noise term in the i' cell as 



" i = W 



' N N 



\j=i 
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(15) 



where A^ is the number of galaxies in the cell. 

For a minimum variance estimator, we wish to choose 
the weights, Wi such that equation (114) is minimised with 
respect to Wi . Denoting the true intrinsic ellipticity variance 
as "true-: this optimal weighting scheme is given by 



= l" 2 ("lue) + cr 2 (A/i)] _1 = (2a? rue + 2M 2 )~ 



(16) 



where " 2 (" 2 rue ) is the uncertainty on the true ellipticity vari- 
ance, and a 2 {Mi) is the uncertainty on the estimated noise 
term. Assuming a 2 rue to be Gaussian distributed, the un- 
certainty on the true ellipticity variance will be " 2 (" 2 rU e) = 
"true- The uncertainty in the noise term can be estimated by 
noting that the noise contribution from randomly orientated 
galaxies will be Mi = e 2 ms /N where e rms is the root mean 
squared random ellipticity of the galaxies. The variance of 
the noise is then given by a 2 (Mt) — (e^)/N 2 — ef ms /N 2 
where e g is the ellipticity of the individual galaxies in cell 
i. Assuming a Gaussian distribution for the galaxy elliptici- 
ties, we can make the approximation, (e^) fa 3(e 2 ) 2 = 3ef m3 
and we have a 2 (Mi) — 2e^ m3 /N 2 = 2M 2 . Combining these 
variances yields the weighting in equation (|16|). Substituting 
this into the expression for a 2 nt yields the minimum variance 
estimator: 



2 _ YsM-^l^true+tM 2 ) 
EP^true+^M 2 )- 1 ' 

The error in this expression is given by 



(17) 



2/ 2 n 
" ("int) — 



i available on tittp: //www. ast.cam.ac.uk/~mike/apmcat/ 



yV = Y.( 2 ^ue + 2M 2 )' 1 . (18) 

Equations (jlTl) and (h8|) form the basis of our analysis. 

4 RESULTS 

4.1 Measurement of the ellipticity variance 

We have used equation ( |l7| ) with atrue = 0.01 to calculate 
the variance of the intrinsic cell ellipticities, af nt over a wide 
range of angular scales. Even so, the results are fairly insen- 
sitive to the value of atrue used. We have measured both 
the bj band variance, (ej(bj)ei(bj)), and the R band vari- 
ance, (ei(R)ei(R)) as well as the cross-correlation between 
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Figure 10. Measured cllipticity variance, af ' . over 436 Supcr- 
COSMOS fields as a function of angular scale, 9 along with the 
weak lensing prediction (heavy dashed line) for a ACDM model 
and a median redshift equal to that of our galaxy sample. Also 
shown are various predictions for the intrinsic signal - see text for 
details. 



the two bands, (ei(bj)ei(R)) . Here, the angled brackets de- 
note the weighted average as in equation (|l7|). Measuring 
the cross-correlation has the advantage that systematic ef- 
fects that are uncorrelated between the two bands are can- 
celled out. The results of all these measurements are shown 
in Figure [id. In the final analysis, we have made use of 
1.86 x 10 6 galaxies for the bj band variance, 1.97 x 10 6 for 
the R band variance and 1.68 x 10 6 for the covariance mea- 
surement. Note that the scale at which the stellar ellipticities 
are smoothed for the PSF anisotropy correction is 15' and so, 
at scales smaller than this, our results may be compromised 
by residual PSF anisotropy distortions. We have, therefore, 
not plotted results below this scale in Figure llOJ. Beyond 
this scale, we are confident that our measurements are not 
dominated by systematics. This assertion is supported by 
the agreement between the two single band variance mea- 
surements and the cross-correlation signal. Further tests for 
systematics are presented in Section 4.2. 

In addition to using equation (J18J) , we have also esti- 
mated the errors on our measurements from 1000 randomi- 
sations of the dataset. In these randomisations, only galaxies 
used in the final analysis were included. The galaxies were 
assigned a position angle taken from a random distribution 
between and 180 degrees while their semi-major and semi- 
minor axes were left unchanged. Equation (n7|) was then 
used to measure the ellipticity variance from the 1000 ran- 
domised datasets. The resulting distributions of 1000 vari- 
ance measurements are centered on zero at all scales, indi- 
cating that our estimator is unbiased. These distributions 
are shown in Figure O, along with the measured signal, for 
four scales: 19, 40 and 56 and 93 arcmins. The errors on the 
signal measured from the real dataset were then calculated 
as the standard deviation of the measurements from the ran- 
domised catalogues. These errors agree quite well with the 
errors calculated using equation (jlq) although those taken 
from the randomisations are, in general, slightly larger. We 
have therefore plotted the la errors from the randomisations 
in Figure IQ. 





o ($=1 8.65 arcmins) 



-5x10 5 5x10 5 10 4 

o 2 (6=39.96 arcmins) 




[jjfj] 


Ml 


i 



-5x10 3 5x10 D 

o (0=55.95 arcmins) 



5x10 

a 2 (S=93.25 arcmins) 



Figure 11. The variance measurements from the 1000 ran- 
domised catalogues (histograms) along with the measured signal 
for four different values of the cell side length, 8. In each case, the 
full line is the bj band variance, the dashed line is the R band 
variance and the dotted line is the covariance measured between 
the two bands. 
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Figure 12. Top panel: The average galaxy ellipticity components, 
(e a ) 2 over the 436 SuperCOSMOS plates for the bj and R bands. 
Bottom panel: The star-galaxy covariance, (e*e;) measured over 
the whole survey for the two bands. 

Our measurements are two orders of magnitude larger 
than the signal expected from weak lensing (see e.g. Jain & 
Seljak, 1997) for a median source redshift of z ~ 0.1 cor- 
responding to the magnitude cut of bj < 20.5 which we 
have used in our analysis. In Figure hd we plot our results 
along with the predicted weak lensing signal for a cluster 
normalised ACDM model, with fi m = 0.3 and $!a = 0.7. 
The discrepancy is exhibited over the entire range of angu- 
lar scale, 9 suggesting that we have not measured extrinsic 
gravitational lensing. 



4.2 Tests for systematics 

To test our results for internal systematics we have esti- 
mated the means of the cell ellipticity components averaged 
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Figure 13. The signal measured (points with error bars) in bj 
(top), R (middle) and the cross-correlation between the two bands 
(bottom) along with the cross correlation of e\ with e 2 for the 
galaxies, (e\e%) (joined crosses). Note the much higher signifi- 
cance of the shear covariance measurement, (e;(bj)ei(i?)} com- 
pared with the shear variance measurements from the individual 
bands alone. 

over all plates. In Figure [12j we show that the mean fields 
are negligible {(e a ) 2 < 10~ ) on all scales. We have also 
estimated the star-galaxy covariance, 



(ee*) = (eie*) + (e 2 ej), 



(19) 



where e* is the stellar ellipticity. This is also shown in Figure 
\L& and is well below our results (|(ee*)| < 10 -5 ) on all 
scales, indicating that, after applying the cut on galaxy size 
(Section 3.2), galaxy ellipticities are no longer correlated 
with stellar values. We have estimated the cross-correlation 
of ei and e 2 for the galaxies, (eie2). This is plotted in Figure 
|l3| along with the measured signal for the two single band 
variance measurements as well as for the shear covariance 
signal. In all cases, the galaxy cross correlation, (eie2) is 
consistently below the measured signal. In particular, for 
the covariance measurement, we have estimated the galaxy 
cross correlation as 

(eie 2 > = <(e 1 (b J )e 2 ( J R) + ei(fl)ea(bj))/2) (20) 

and this is consistent with zero at the 10 -9 level on all scales, 
indicating that (ei(bj)ei(R)) is free of significant residual 
systematics and is a robust estimate of the galaxy alignment 
signal. 

As a test for systematics introduced by the measuring 
process we have pushed the APM data (Section 3.5) through 
our analysis for one plate. In Figure [14 we plot the variance 
as measured by the two machines for this field. Again we 
find a very good agreement between both catalogues, sug- 
gesting that any systematic effects in our analysis are small, 
°"sys < 10 -5 , and below our measurement. This agreement 
between the catalogues is demonstrated in Figure [13 where 
the difference in variance measurements from the two sur- 
veys is plotted. This difference is consistent with zero on all 
scales. 

As noted in Section 3.1, the SuperCOSMOS machine 
measures unweighted second moments of all objects detected 
within an isophotal threshold, I t h,- The nominal value of this 
threshold in the SuperCOSMOS Sky Survey is 2.3u above 



e(arcmins) 
Figure 14. Comparison of the SuperCOSMOS variance of ellip- 
ticities and the APM variance for one plate (UKST field 78). Note 
the errors are much larger than in the final analysis due to the 
much smaller number of galaxies measured. The SuperCOSMOS 
points have been slightly laterally displaced for clarity. 




0(arcmins) 
Figure 15. The difference in the ellipticity variance as measured 
from the SuperCOSMOS and APM machines (Act 2 = o- 2 pm - 
a scos ) f° r * ne same field as shown in figure h 4. The value of Act 2 
is consistent with zero on all scales indicating that systematic 
errors introduced by the plate scanning procedure are small. 

the sky background, I s k y - To further establish the reliability 
of our ellipticity measurements, we have varied I t h. about 
it's nominal value for one field and compared the resulting 
galaxy ellipticity measurements with those from the main 
analysis for this field. The results of this test are shown 
in Figure hq showing that the determination of the galaxy 
ellipticities is largely insensitive to the isophotal threshold 
used to measure the object parameters. 



4.3 Intrinsic alignments 

Recent numerical and analytic studies (Catelan et al., 2001; 
Crittenden et al., 2001; Croft & Metzler, 2001; Heavens et 
al., 2000) have put constraints on the shear signal expected 
from intrinsic galaxy alignments. In these studies the galaxy 
shape is estimated either from the perpendicular to the halo 
angular momentum vector for spirals, or from the shape of 
the halo for ellipticals. 
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Figure 16. Individual galaxy cllipticities as measured from 
UKST field 78 with different isophotal thresholds applied. In each 
of the four panels, the horizontal axis measurements are from the 
Sky Survey data with an isophotal threshold value of 2.3ct above 
the sky background. The vertical axes measurements have been 
made using isophotal threshold values of (a) 1.5ct, (b) 2. Oct, (c) 
3.5ct and (d) 4.5ct above the sky background. 



In Figure hoi we have plotted our results along with 
the alignment effect estimated from each group. Also shown 
are results from an extension to the numerical simulations 
work of Heavens et. al., 2000 for the spiral galaxy model at 
Z — 0.1 (A. Heavens, private communication). The spread 
in models can be accounted for by the choice of mechanism 
and galaxy type, assumptions about alignments between ha- 
los and galaxies, correlations between tidal and inertial fields 
and disc thicknesses. We have used the redshift scaling z~ 2n , 
suggested by Crittenden et al., where n is the slope of the 
matter correlation function, to scale the results from differ- 
ent redshifts. We have also assumed a 2 (9) = 2/(2 — n)C{6), 
where C(9) is the ellipticity correlation function. 

Generalising the arguments by Crittenden et al. (2001), 
at large scales the ellipticity variance should scale as 



a 2 (6) 



Az 



- 2n |l + (cV0 o ) 2 | 



(21) 



One would hope that with future, more precise observations 
of intrinsic alignments, it may become possible to constrain 
the slope of the matter correlation function, and thus shed 
light on the clustering properties of dark matter on these 
scales (< lO/i^Mpc). 



5 CONCLUSION 

In this paper we have presented a measurement of the intrin- 
sic alignment effect of galaxies on scales from a few arcmin- 
utes to 100 arcminutes. Using 2 million galaxy ellipticities 
measured from the digitized SuperCOSMOS Sky Survey, to 
a depth of bj = 20.5 and median redshift z — 0.1 and cover- 
ing 436 plates or 10,000 sq. deg. in two passbands, we have 
corrected the data for distortions due to PSF anisotropy and 
seeing, and applied a new minimum variance estimator to 
the data. 

After applying these corrections to the dataset and 



excluding the over-corrected small galaxies, we find our 
measurements are internally consistent, with good agree- 
ment in regions where the plates overlap, with effectively 
zero star-galaxy ellipticity correlations and zero ei-e2 cross- 
correlations. The ellipticity variance measurements from the 
bj and R passbands agree well with each other and with the 
cross-correlation between the two bands — further confirma- 
tion that systematic effects are small. We have also demon- 
strated external consistency with the APM sky catalogue 
data over a restricted region of sky, indicating that we are 
not contaminated by measurement systematics. 

Our resulting estimates of the ellipticity variance over 
a wide range of scales are two orders of magnitude higher 
than that expected from gravitational lensing by large-scale 
structure, but roughly in line with those predicted from in- 
trinsic alignments in the gravitational instability scenario, 
although the predictions for intrinsic alignments are still un- 
certain. For instance, it is not clear if the effect is dominated 
by tidal spin or shape alignments. This agreement suggests 
that we are not contaminated by internal systematics. 

Our results imply that other shallow surveys, such as 
the Sloan Digital Sky Survey, should measure roughly the 
same contribution to the total variance from intrinsic align- 
ments and gravitational lensing shear. Since the intrinsic 
alignment signal is expected to scale as z~ 2 (Crittenden et 
al., 2001) and the lensing signal as z (Jain & Seljak, 1997), 
we expect the ratio of intrinsic alignment to gravitational 
shear to scale as 



^- ^10 2 



(£)' 



(22) 



This ratio is unity at around z ~ 0.37. Higher redshift sur- 
veys, such as the VISTA Gravitational Lensing Survey (Tay- 
lor, 2001) with z k 1, will measure gravitational lensing 
shear. 

While one would hope that the correlations between el- 
lipticity, spin and the shear field would allow one to measure 
the amplitude of the dark matter density field, the spin is 
also determined by the inertial tensor of the dark matter 
halo. At present, the relation between shape and the local 
shear field of haloes, and the relationship between galaxy 
ellipticity and halo shape, introduces a large uncertainty in 
our understanding of intrinsic ellipticity alignments. With 
present and future observations and theory we can hope to 
resolve these issues. 
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